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Abstract 

Previous  attempts  at  modeling  sound  propagation  through  a  forest  have  largely  neglected  the  effects 
of  a  sound  speed  profile.  This  paper  presents  a  PE-based  sound  propagation  model  that  includes 
forest  effects.  In  addition  to  a  simplified,  but  realistic,  sound  speed  profile,  the  model  includes  ground 
impedance  effects,  bulk  attenuation  due  to  mulitple  scattering  by  the  tree  trunks  and  canopy,  and  the 
usual  spherical  spreading  and  atmospheric  absorption.  The  model  is  compared  to  acoustic  measurements 
taken  in  July  2002  at  Lone  Star  Army  Ammunition  Plant  in  Texarkana,  TX.  In  the  measurement, 
Composition  C-4  explosions  were  used  as  the  acoustic  source.  Acoustic  measurements  were  make  at 
174  m  and  1400  m  from  the  source  in  the  forest,  and  174  m  from  the  source  in  an  open  field.  Measured 
tree  trunk  diameters,  number  of  trees  per  acre,  tree  height,  canopy  height,  and  ground  impedance  are 
used  as  parameters  in  the  PE  model.  A  predicted  C-4  spectrum  is  used  to  weight  the  starting  field  in  the 
PE.  Upwind  and  downwind  cases  are  examined.  The  forest  PE  model  is  compared,  using  appropriate 
parameters,  to  both  the  forest  and  the  open  field  measurements.  Finally,  the  forest  PE  model  is  compared 
to  the  open  field  PE  model  for  long  distances,  as  a  way  to  investigate  differences  in  propagation  through 
these  distinctive  settings. 


Part  I 

Introduction 


The  forest  is  a  highly  complex  acoustic  medium.  In  attempting  to  understand  it,  one  must  consider  the 
vertical  sound  speed  profile,  the  ground  impedance,  and  scattering  from  the  trees.  Previous  attempts 
at  modeling  the  forest  have  typically  considered  only  two  of  the  three  major  components,  usually  ground 
impedance  and  tree  scattering,  although  more  recent  papers  included  only  the  vertical  sound  speed  profile 
(see,  for  example,  Swearingen  [12],  Albert  [1],  Tunick  [13],  and  Heimann  [9]).  In  this  paper  we  present  a 
model  that  incorporates  all  three  factors.  The  model  is  based  on  the  GFPE  algorithm.  The  model  is  used 
in  two  different  ways,  first  to  compare  forested  and  open  field  propagation  to  experimental  data,  and  then 
to  compare  forested  to  open  field  propagation.  These  comparisons  are  performed  at  174  m  and  1400  m, 
because  these  distances  correspond  to  experimental  receiver  locations. 
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Part  II 

The  Forest  Model 

As  stated  in  the  introduction,  the  model  is  based  on  the  GFPE.  Since  that  model  is  well  described  in 
the  literature,  the  derivation  will  not  be  included  here.  Instead,  the  reader  is  referred  to  Gilbert  [8]  and 
Salomons  [11],  two  excellent  papers  describing  the  GFPE.  The  following  sections  will  describe  the  geometry, 
vertical  sound  speed  profile,  ground  impedance,  trunk  and  canopy  scattering,  and  the  blast  spectrum. 


1  Geometry 

The  geometry  of  this  model  is  based  on  measurements  taken  in  the  field  test.  This  dictates  that  the  forest 
height  be  15  m,  and  the  canopy  begins  at  10  m  above  the  ground.  While  the  trunk  and  canopy  sizes  and 
spacings  are  variable  within  the  model  itself,  they  are  set  to  match  the  measurements  taken  in  the  field.  The 
source  height  was  2  m  above  the  ground,  and,  for  the  spectra,  the  receiver  height  is  1.2  m  above  the  ground. 
Trunk  cylinders  are  randomly  spaced,  but  all  perpendicular  to  the  ground,  and  extend  to  the  top  of  the 
canopy.  The  canopy  is  also  made  up  of  cylinders,  which  are  oriented  parallel  to  the  trunks.  They  extend 
from  the  bottom  of  the  canopy  to  the  top  of  the  canopy.  Only  two  estimated  branch  sizes  are  considered. 
Figure  1  is  a  graphical  interpretation  of  the  model. 


Figure  1 :  Geometry  of  the  model 
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2  Atmospheric  Variation  with  Height 

The  atmosphere  near  the  ground  surface  exchanges  heat  with  the  surface  and  loses  speed  due  to  viscous  drag. 
This  drag  causes  shear  in  the  wind  profile  and  is  the  major  generator  of  mechanical  turbulence.  Depending 
on  the  surface  heating,  the  wind  speed  variation  with  height  above  the  boundary  can  take  one  of  several 
forms.  An  excellent  recent  summary  of  the  situation  has  been  published  by  Wilson  [15]. 

The  presence  of  the  forest  canopy  serves  as  a  target  for  solar  heating,  causing  the  ground  to  be  partly 
shaded.  In  the  absence  of  solar  heating,  radiative  cooling  within  the  canopy  is  reduced  as  well,  because  the 
upper  canopy  is  at  a  much  higher  temperature  than  the  open  night  sky.  The  presence  of  trees  appears  to 
the  wind  as  additional  roughness  elements,  increasing  the  drag.  Because  leaf  areas,  heat  capacities,  and 
transpiration  of  the  various  surfaces  differ  widely  between  tree  types,  the  detailed  microclimate  within  the 
forest  is  complex.  Tunick[13]  has  developed  a  two-dimensional  model  for  this  problem. 

In  order  to  understand  refraction  of  traffic  noise  beyond  a  grove  of  trees,  Heimann[9]  implemented  a  3-D 
computational  fluid  dynamical  model  for  the  flow  past  a  random  assembly  of  truncated  vertical  right  circular 
cylinders.  The  supposed  incident  flow  field  consisted  of  a  wind  profile  with  logarithmic  variation  with  height. 
After  passing  through  the  grove  of  trees,  the  mean  flow  at  mid-height  of  the  canopy  was  reduced  and  nearly 
constant  with  height.  Near  the  top  of  the  canopy  was  a  strong  wind  shear,  with  wind  increasing  a  bit  more 
than  the  incident  wind  at  that  height.  Above  the  canopy,  the  difference  between  forest  wind  and  open  field 
wind  grew  smaller. 

Because  the  incident  open  field  wind  and  the  forest  wind  profiles  result  from  the  same  impressed  wind 
field,  they  offer  a  foundation  for  comparative  assessment  of  wind  refraction.  We  chose  to  use  the  wind  fields 
as  computed  by  Heimann[9]  as  a  starting  point  for  examining  the  basic  effect  of  wind  refraction  of  sound  by 
forest  and  open  field.  We  scaled  the  wind  profiles  by  height  to  match  the  canopy  height  in  our  experimental 
measurement,  and  extended  the  published  profiles  to  smoothly  converge  to  same  wind  speed  at  40  m  height, 
and  to  conserve  total  horizontal  flow.  Note  that  just  above  the  canopy  top,  the  comparative  increase  of  wind 
speed  above  the  open  field  value  is  required  to  balance  the  surplus  of  incident  wind  from  the  open  field  at 
the  canopy  height  and  below.  The  scaled  wind  profiles  are  shown  in  Figure  2. 

Of  course,  left  unaccounted  for  by  this  scaling  were  adjustments  for  the  tree  stand  densities  and  the  leaf 
area  indices.  We  have  also  omitted  the  temperature  dependence  of  the  problem  by  setting  the  temperature 
to  a  constant.  Nevertheless,  the  strengths  of  the  gradients  are  appropriate  for  some  value  of  imposed  flow. 
Additionally,  the  larger  wind  profile  gradients  coincide  with  the  ground  surface  and  canopy  top,  as  would 
temperature  gradients  resulting  from  surface  radiative  heat  transfer. 
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Hind  Velocity  Profile 


Figure  2:  Model  wind  profiles  for  forest  and  open  field 


3  Ground  Impedance 


Ground  impedance  under  the  tree  canopy  is  typically  affected  most  by  the  detritous  composing  the  surface 
layer.  The  ground  impedance  was  not  directly  measured.  Instead,  the  measured  receiver  spectrum  at  short 
range  (<  80  m  from  the  source)  was  compared  to  a  forward  prediction  using  a  van  der  Pol-type  solution  for 
a  homogeneous  half-space  above  porous  ground.  The  ground  impedance  was  modeled  with  a  simple  function 
developed  by  Attenborough  [3],  providing  a  reasonable  match  to  the  receiver  spectrum.  The  ground  "dip" 
in  the  spectrum  was  fairly  easy  to  obtain  by  adjusting  effective  flow  resistivity  until  good  fit  to  the  receiver 
spectrum  was  achieved. 

Within  the  trees,  the  porous  surface  was  made  up  of  a  layer  of  pine  straw  (3  cm),  above  a  very  thin 
(<  0.3  cm)  layer  of  highly-decomposed  pine  straw.  Below  this  was  a  fairly-consolidated  sand  layer  (with  no 
evidence  of  decomposed  straw),  which  was  modeled  as  a  hard-backing. 

In  the  open  field,  a  mixture  of  thick  grasses  (lm  high)  and  blackberry  bushes  (lm  high),  overlaid  a 
spongy  surface  atop  clay.  A  few  points  on  the  surface  were  saturated  with  water;  others  were  dry.  The 
surface  was  sloped  approximately  1  to  2  degrees  transverse  to  the  direction  of  propagation. 

The  characteristic  impedance  of  the  porous  surface  layer  is  calculated  from  Zc  —  upb/kb,  as  a  fluid  having 
complex-valued  bulk  density  and  wavenumber, 
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where  T  is  tortuousity,  O  is  volume  porosity  (m3/m3),  u  is  angular  frequency  (rad/s),  p  is  air  density 
(kg/ m3),  7  is  the  ratio  of  specific  heat  capacities,  and  Npr  is  the  Prandtl  number. 

This  ground  surface  corresponds  to  Attenborough’s  three  parameter  homogeneous  (H3A)  model.  Finally, 
using  the  H3A  model  for  a  thin  upper  layer  where  porosity  decreases  with  depth,  and  with  hard  backing, 
leads  to 

Z/ pc  =  (Zcj pc)  +  i/  Q'yupQ.d) 
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where  c  is  the  speed  of  sound  in  air  ( m/  s),  and  d  is  the  depth  of  the  surface  layer  ( m). 

The  parameters  used  for  evaluating  ground  impedance  with  the  H3A  model  were  p  =  1.2,  c  =  350, 
7  =  1.4,  and  Npr  =  0.724.  For  pine  forest:  cre  =  7.5 kPa  •  s  •  m“2,  d  =  0.03  m,  0  =  0.4,  and  T  =  2.0.  For 
open  field,  ae  =  57kPa  •  s  *  m-2,  d  =  0.053m,  Cl  =  0.4,  and  T  —  2.0. 


4  Trunk  and  Canopy  Scattering 

The  trunk  and  canopy  scattering  are  both  based  on  Twersky’s  multiple  scattering  model  (see,  for  example 
Twersky  [14]  or  Embleton  [5]).  In  both  the  trunk  and  canopy  cases,  the  approximation  of  infinite  cylinders 
is  used,  scaled  appropriately  for  size  and  number  density  of  scatterers.  The  number  and  size  of  tree  trunks 
are  taken  from  the  measurements.  The  number  density  of  trees  is  0.0124  m”1  and  the  average  tree  radius  is 
0.0925  m.  The  canopy  is  approximated  by  large  and  small  branches.  These  are  scaled  in  size  and  number 
from  the  tree  trunk  parameters.  The  large  branches  are  assumed  to  be  1/3  of  the  trunk  radius  and  six 
times  the  number  of  trunks.  The  small  branches  are  assumed  to  be  1/8  of  the  trunk  radius  and  24  times 
the  number  of  trunks. 

The  result  of  Twersky’s  multiple  scattering  model  for  randomly  placed  parallel  infinite  cylinders  is  a 
complex  effective  bulk  propagation  wavenumber,  ksi  where 


ikl  -  4 iNg  +  (; g '2 


-g*)(2N/k0)2 


and  k0  —  uj/c0  is  the  unobtruded  wavenumber,  N  is  the  number  of  trees  per  meter,  and  g  and  g'  are  the 
forward  and  backward  scattering  coefficients,  respectively,  defined  as 
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and  Zw  is  the  surface  impedance  of  the  tree.  In  this  case,  Zw  is  nearly  rigid,  following  Price[10].  Plots  of 
the  real  and  imaginary  parts  of  ks  are  included  in  Figure  3. 

The  scattering  is  implemented  in  the  GFPE  by  adjusting  the  wavenumber  in  the  appropriate  regions. 
For  the  trunk  scattering,  the  wavenumber  is  adjusted  from  the  ground  to  the  tree  top  height,  or  0  m  to 
15  m.  The  canopy  begins  at  10.6  m.  To  incorporate  the  canopy  scattering,  the  imaginary  parts  of  the  bulk 
scattering  wavenumbers  due  to  the  large  and  small  branches  are  added  to  the  trunk  scattering  wavenumber. 
This  is  possible  because  the  real  part  of  the  wavenumbers  are  constant  within  0.02%,  and  are  therefore 
considered  negligible  differences  (see  Figure  4). 

While  the  scattering  effect  is  small  over  short  distances  and  for  low  frequencies,  it  increases  significantly 
with  both  frequency  and  distance.  Figure  5  shows  the  attenuation  per  unit  thickness  due  to  scattering  in 
units  of  dB/lOOm.  The  plot  shows  that  attenuation  increases  with  frequency,  as  expected,  and  that  it  has 
a  greater  overall  effect  than  atmospheric  absorption. 

Atree  =  20  •  Im(fcs)  •  distf  log10(e)  (dB)  (1) 
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Figure  3:  Real  and  imaginary  parts  of  ks  plotted  as  a  function  of  frequency 
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Figure  4:  Relative  differences  between  scattering  wavenumbers  and  unobtruded  wavenumber 
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Attenuation  rate  due  to  scattering 
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Figure  5:  Attenuation  rate  due  to  scattering,  in  dB/100  m 

5  Blast  Wave  Spectrum 

Explosions  of  Composition  C-4  were  employed  as  the  sound  source.  The  use  of  C-4  provided  a  high-pressure, 
broadband,  compact,  isotropic  source  with  fairly  repeatable  event-to-event  sound  energy.  The  charges  were 
formed  into  spheres,  each  wrapped  with  a  thin  plastic  bag  and  suspended  by  a  rope  so  that  their  centers  were 
approximately  2  m  above  the  surface.  The  sound  wave  produced  by  an  explosion  in  air  has  been  modeled 
extensively.  It  is  convenient  to  approximate  the  sound  pressure-time  history  with  an  expression  due  to 
Friedlander  [7], 

pt  =  p°(l-  ^)  exp  (-t/td)  &t  (2) 

where  p°  is  the  peak  positive  overpressure  ( Pa),  td  is  the  positive  phase  duration  ( s),  and  ©  is  the  Heaviside 
step  operator. 

Although  the  waveform  fails  to  indicate  the  finite  rise  time  of  the  shock,  or  to  show  pulse  broadening  or 
rounding  that  result  from  attenuation  and  temporal  dispersion,  it  offers  a  simple  starting  point  for  spectral 
analysis.  Its  Fourier  transform  is, 

pf  =  f  pte~lu>tdt  =  pQtd  (—iu>td)  /  (1  +  iu)td)2  (3) 

J — oo 

where  u  =  2tt/  is  the  angular  frequency,  and  it  has  the  corresponding  sound  exposure  spectrum, 

Ef  =  2  (p°td)2  u2t2d/ (1  +  u2t2d)2 .  (4) 

This  function  has  its  maximum  value  at  a;  =  l/td. 

For  the  (admittedly  unrealistic)  homogeneous  and  unbounded  field  at  1km  distance  from  a  1kg  charge 
of  Composition  C-4,  p\ kg  =  79 Pa  and  td, ikg  =  0.0089s.  For  other  charge  sizes,  it  generally  holds  that  the 

pressure  and  duration  increase  according  to  the  cube  root  of  charge  weight,  p^/Pikg  “  (W,/ 1  kg)1/3  and 

td,x/td, ikg  =  (wx/l  kg)1/3,  although  more  accurate  expressions  appear  in  ANSI  S2.20[2]  and  Ford  et  al.  [6]. 

Figure  6  shows  the  spectrum  of  a  C-4  explosion  at  1  km  in  an  unbounded  field,  evaluated  for  the  two 
charge  sizes  employed  in  the  measurement,  0.57  kg  and  2.27  kg.  Both  narrow  band  (1Hz  bandwidth)  and 
1/3-octave  spectra  are  shown.  Also  shown  axe  curves  calculated  by  adding  101g(/c/lHz)  —  6.36  dB  to  the 
narrow  band  sound  exposure  levels.  These  latter  curves  are  shown  to  enable  rough  comparison  to  1/3-octave 
band  spectra. 
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Spectral  Sound  Exposure  of  C-4  Explosion  at  1  kn:  Unbounded  Field 


Figure  6:  Idealized  sound  exposure  spectrum  of  Composition  C-4  explosion  at  1km  distance  in  an  unbounded 
homogeneous  atmosphere 


6  The  Forest  GFPE 


The  previous  sections  all  described  the  components  of  the  forest  GFPE  model.  This  section  contains  example 
plots  of  the  resulting  model.  The  model  was  calculated  for  both  upwind  and  downwind  propagation  over 
a  range  of  frequencies  from  1Hz  to  2200  Hz.  The  figures  below  are  snapshots  taken  at  40  Hz,  400  Hz,  and 
1000  Hz.  In  all  of  the  figures,  the  source  is  located  at  2  m  above  the  ground,  and  is  assumed  to  be  a  unit 
source  with  an  analytic  starting  field,  as  shown  below.  The  form  of  the  starting  field  follows  that  in  Cooper[4] 


<Kra,z)  =  { 


+  eJw[R* + (1  - 

0, 


0  <  z  <  z3 

Z8  <  Z  <  Zattn 
otherwise 


where  ra  =  0  +  dr,  i?i  =  ^(rsf  +  (z-  zs)2,  R2  =  \/(r,)2  +  (z  +  zs)2,  Rp  is  the  plane  wave  reflection 
coefficient,  given  by 


Rp  =  (sin  (j>  -  Pg)/ (sin  <f>  +  Pg ) 

where  P  is  the  normalized  acoustical  admittance  of  the  ground  surface,  pg  =  l/Zg,  and 


Fi  =  1  +  j\fnh.e  erfc(-jA) 

A  =  yf\ ikR2  yjl  +  Pg  sin  0  -  (1  -  cos  <j>  (5) 


where  <j>  is  the  acute  angle  subtended  by  R2  and  the  ground,  and  Fx  represents  the  first  term  in  the  asymptotic 
expansion  solution  of  the  spherical  wave  reflection  coefficient. 

In  all  of  the  figures,  notice  that  the  forest  appears  to  be  "trapping"  energy  below  the  canopy.  This  effect 
is  much  more  pronounced  in  the  downwind  cases.  Shown  below  are  the  following:  40  Hz  in  Figure  7,  400  Hz 
in  Figure  8,  and  1000  Hz  in  Figure  9. 

We  will  next  examine  the  relative  contributions  of  each  component  of  the  model  by  starting  with  only 
the  atmospheric  profile  present,  and  then  adding  additional  components.  This  has  been  done  in  Figures  10 
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and  11,  both  upwind  and  downwind  cases.  Looking  at  the  174  m  location,  the  ground  effect  is  large,  and  the 
scattering  components  are  only  roughly  a  2  dB  contribution  at  1000  Hz.  At  the  1400  m  location,  the  ground 
effect  is  again  large,  but  the  scattering  is  now  much  more  significant,  especially  above  around  400  Hz. 


Upwind,  40  Hz 


Range  m 


Figure  7:  Forest  GFPE  at  40  Hz,  upwind  and  downwind  cases. 
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Height  m  Height  m  Height  m  Height  m 


figure  8:  Forest  GFPE  at  400  Hz,  upwind  and  downwind  cases 


Upwind,  1000  Hz 


Range  m 


Figure  9:  Forest  GFPE  at  1000  Hz,  upwind  and  downwind  cases 


Range  m 


Upwind  Forest,  174  m 


Frequency  (Hz) 
Downwind  Forest,  174  m 


Frequency  (Hz) 


Figure  10:  Buildup  of  model  components,  174  m  location.  In  the  legend,  A= Atmospheric  Profile,  G=Ground 
Effect,  T— Trunk  Scattering,  C=Canopy  Scattering. 


Upwind  Forest,  1400  m 


Figure  11:  Buildup  of  model  components,  1400  m  location.  In  the  legend,  A= Atmospheric  Profile, 
G=Ground  Effect,  T=Trunk  Scattering,  C^Canopy  Scattering. 
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Part  III 

Analysis 

The  forest  GFPE  has  been  analyzed  two  different  ways:  by  comparisons  to  data  and  comparisons  between 
forest  and  open  field.  These  comparisons  are  performed  in  both  short  (174  m)  and  long  (1400  m)  range 
regimes. 


7  Comparisons  to  Data 

Comparing  theoretical  predictions  to  field  data  can  lead  to  interesting  observations  about  both  the  prediction 
and  the  field  test.  In  this  case,  data  from  a  open  field-forest  experiment  was  available.  In  the  experiment, 
two  sizes  of  Composition  C-4  were  detonated  at  a  height  of  2  m  above  the  ground,  from  four  different 
locations.  There  were  two  source  locations  in  the  open  field  and  just  along  the  edge  of  the  forest.  The  only 
data  used  here  was  taken  in  an  open  field  at  a  distance  of  174  m  from  one  open  field  source,  in  the  forest 
174  m  from  one  forest  source,  and  in  the  forest  1400  m  from  one  forest  source.  In  all  cases,  the  measured 
data  is  presented  in  1/3  octave  band  SEL.  The  prediction  spectra  are  weighted  by  the  appropriate  blast 
source  spectrum. 

Figure  12  shows  the  comparison  between  the  forest  GFPE  calculations  and  the  measured  forest  propa¬ 
gation  data  at  174  m  from  the  source.  At  short  range,  the  results  are  heavily  dependent  on  assumptions  we 
have  made  about  the  source.  It  is  possible  that  the  source  spectrum  used  is  not  the  most  appropriate  one. 
The  measured  data  is  generally  higher  in  level  than  the  predictions.  The  ground  effect  in  the  model  is  more 
pronounced  than  in  the  measurements.  However,  the  general  shapes  match  reasonably  well.  The  ground 
dip  and  maximum  value  both  appear  at  the  proper  frequencies.  Further  improvements  are  anticipated  with 
adjustments  to  parameters,  such  as  ground  impedance. 

Figure  13  shows  the  comparison  between  the  field  GFPE  calculations  and  the  measured  open  field  prop¬ 
agation  data.  Again,  at  short  range,  the  model  results  are  heavily  dependent  on  the  source  model.  The 
apparent  ground  dip  in  the  downwind  prediction  is  not  seen  in  the  measurement. 

Figure  14  shows  the  1.251b  C-4  case  for  a  fully  forested  propagation  path  at  1400  m.  After  more  careful 
analysis  of  the  measured  data,  we  determined  that  the  spectrum  in  this  comparison  is  likely  to  be  anomalously 
low  in  level.  This  conclusion  was  derived  from  the  fact  that  most  of  the  long  range  data  was  uncalibrated  or 
electronically  clipped  and  could  not  be  used  due  to  equipment  errors.  The  measured  spectrum  shown  is  from 
an  un-clipped  signal,  and  is  therefore  likely  to  be  a  lower  level  than  the  others.  However,  it  is  important  to 
notice  that  the  general  shape  of  the  spectrum  matches  reasonably  well  with  the  predicted  spectra. 
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Forest,  174m,  1.251b  C-4 


Frequency  (Hz) 


Forest,  174m,  5tb  C-4 


Frequency  (Hz) 


Figure  12:  Comparison  of  forest  measurements  to  GFPE  model,  174  m. 


Reid,  174m,  1.251b  C-4 


Frequency  (Hz) 


Reid,  174m,  51b  C-4 


Frequency  (Hz) 


Figure  13:  Comparison  of  open  field  measurements  with  GFPE  model,  174  m. 
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Forest,  1400m,  1.251b  C-4 


Frequency  (Hz) 


Figure  14:  Comparison  of  forest  measurements  to  GFPE  model,  1400  m. 


8  Forest  vs.  Field  Comparisons 

The  original  intent  of  the  development  of  this  model  was  to  determine  whether  the  presence  of  a  forest 
around  a  noise  source  actually  helped  to  mitigate  noise.  One  way  to  get  at  this  information  from  the  model 
is  to  compare  open  field  predictions  to  forest  predictions.  Figure  15  shows  the  difference  in  level  between 
open  field  and  forest  predictions.  If  the  resulting  difference  is  positive,  the  open  field  is  a  better  propagation 
medium  over  that  frequency  range.  If  the  resulting  difference  is  negative,  the  forest  is  a  better  propagation 
medium  over  that  frequency  range.  It  is  important  to  note  that  the  model  does  not  include  atmospheric 
turbulence.  Omitting  atmospheric  turbulence  causes  the  open  field  upwind  cases  to  develop  a  strong  shadow 
zone  that  would  be  somewhat  lessened  in  a  real  setting. 

From  the  comparison  plots  in  Figure  15,  several  conclusions  can  be  drawn.  In  the  short  range  (174  m) 
downwind  case,  the  forest  is  beneficial  above  200  Hz,  by  an  average  of  8  dB.  Between  45  and  100  Hz,  there^ 
is  a  slight  benefit  (1-2  dB)  from  the  forest.  Below  45  Hz,  the  forest  is  a  better  propagating  medium  than 
the  open  field.  In  the  short  range  (174  m)  upwind  case,  with  the  exceptions  of  20  Hz  and  300  Hz,  the  open 
field  is  more  beneficial  in  terms  of  noise  mitigation.  In  the  long  range  (1400  m)  downwind  case,  the  forest 
is  beneficial  for  mitigation  above  200  Hz,  on  the  order  of  8-10  dB  on  average,  and  anywhere  from  0-20  dB. 
Below  200  Hz  at  long  range  downwind,  the  open  field  causes  3-11  dB  more  attenuation  than  the  forest.  In 
the  upwind  case,  the  forest  is  more  beneficial  for  mitigation  below  35  Hz,  significantly  worse  between  35  Hz 
and  150  Hz,  and  then  nominally  more  attenuating  than  the  open  field  above  150  Hz. 
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Downwind  Open  -  Downwind  Forest,  174m 


Upwind  Open  -  Upwind  Forest,  174m 


Downwind  Open  -  Downwind  Forest,  1400m 


Upwind  Open  -  Upwind  Forest,  1400m 


Figure  15:  Relative  attenuation  of  forest  vs.  open  field. 


Part  IV 

Conclusions  and  Future  Work 

Overall,  the  forest  GFPE  matches  data  taken  in  the  forest  and  the  open  field  in  shape  and  trend,  although 
not  in  magnitude.  Although  the  levels  are  not  exactly  predicted,  the  general  shapes  of  the  spectra  are 
close,  especially  in  the  long  range.  Possible  sources  of  error  include  the  choice  of  the  Friedlander  spectrum 
to  weight  the  predictions,  the  values  used  for  ground  impedance,  and  the  atmospheric  profile  used.  More 
analysis  is  needed,  including  an  assessment  of  the  model’s  sensitivity  to  various  parameters. 

Because  only  one  case  has  been  examined  here,  it  is  difficult  to  draw  general  conclusions  about  the 
effectiveness  of  a  forest  for  noise  mitigation.  However,  for  the  set-up  presented  here,  the  upwind  open  field 
case  produces  the  lowest  levels  across  almost  all  frequencies  tested,  and  at  both  short  range  (174m),  and  at 
long  range  (1400  m),  with  the  primary  exception  located  near  300  Hz,  which  corresponds  to  the  ground  dip 
in  the  forest.  Again,  it  is  important  to  point  out  that  turbulence  has  been  neglected  in  all  cases.  When 
comparing  downwind  cases,  the  open  field  is  slightly  more  attenuating  than  the  forest  for  frequencies  below 
50  Hz  in  the  short  range,  with  the  forest  providing  more  attenuation  above  50  Hz.  In  the  long  range,  the 
forest  provides  more  attenuation  above  roughly  200  Hz. 

In  the  future,  we  plan  to  improve  the  atmospheric  profile,  perform  sensitivity  testing  on  atmospheric, 
ground  impedance  and  scattering  parameters.  An  important  remaining  challenge  is  to  incorporate  the 
incoherent  field  effects  composing  the  backscatter  and  reverberation  signals. 
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